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Abstract. Coherent diffraction imaging (CDI) on Bragg reflections is a promising 
technique for the study of three-dimensional (3D) composition and strain fields in nano- 
structures, which can be recovered directly from the coherent diffraction data recorded 
on single objects. In this article we report results obtained for single homogeneous 
and heterogeneous nanowires with a diameter smaller than 100 nm, for which we used 
CDI to retrieve information about deformation and faults existing in these wires. The 
article also discusses the influence of stacking faults, which can create artefacts during 
the reconstruction of the nanowire shape and deformation. 
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1. Introduction 

Coherent X-ray Diffraction Imaging (CXDI) has been developped for the past 10 years 
[U I2J El H] : thanks to the development of focusing optics [5] such as compound refractive 
lenses (CRL) [6], Kirckpatrick-Baez mirrors (KB) [7] and Fresnel-zone plates (FZP) 
[9], it is now possible to collect the scattering of single objects with a size smaller 
than 100 nm [10]. 

One goal of CXDI is to allow for the 3D reconstruction of single, non-crystalline 
objects such as biomolecules [TT] or amorphous materials |12J : in this case the small- 
angle scattering is measured, optionally in 3D using a tomographic approach [13], and 
the electronic density of the sample can be reconstructed using an inverse Fourier 
transform combined with phase retrieval algorithms [HI [151 EE fTT] . 

In the case of crystalline materials, it is also possible to measure the scattered 
intensity around Bragg peaks [21 [31 [HI [19] , a method which we will refer to as Coherent 
Bragg Imaging (CBI). CBI allows to recover the shape (the electronic density) of the 
objects, but this is not obviously the main interest of the method, as electron microscopes 
can much more easily give access to the shape of any free-standing nano-object, with a 
better resolution. The interest of CBI comes from the fact that the scattered intensity 
around a Bragg reflection is sensitive to the deformation of the crystal, at a resolution 
smaller than the d-spacing of the considered reflection. Moreover it is possible to study 
buried objects in their normal environment (e.g. beneath a protecting shell) without 
sample preparation that could change their strain state. So far only a few examples of 
deformation reconstruction in nano-objects have been published using CBI [31 [201 l2Tj . 
generally with small deformations that were not due to a controlled epitaxial process. 

Nanowires have recently been the focus of CBI experiments [211 221 E31 EH [25] , but 
so far no CBI on heterogenous nanostructures have been presented with a successful 
reconstruction of the shape and deformation field, although a few measurements on 
SiGe islands have been reported [261 [271 [28] . A possible explanation for the relative lack 
of publication in the field is the existence of structural faults: dislocations (inherent 
to heterogeneous structures grown by epitaxy) and stacking faults (related to the low 
stacking fault energy [29] often found for nanowires). 

In this article we will discuss the importance of structural faults in coherent scatter- 
ing experiments, and how they affect our ability to study heterogeneous nanostructures 
using CBI, in the particular case of nanowires. 

This paper is organized as follows: in section [2] we give a general presentation of 
coherent Bragg imaging in the case of strained nano-objects. Section [3] presents CBI re- 
sults on homogeneous nanowires as well as perfect (simulated) heterogeneous nanowires. 
Section H] details how the presence of stacking faults affects CBI, and how it is possible 
to obtain quantitative results for strain mapping and a statistical analysis of the faults. 

The following crystallographic conventions will be used: in reciprocal space the 
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scattering vector k will be given either in reciprocal lattice units (r.l.u.) (hkl), or using 
its norm k = 2s ™^ . Displacement values (e.g. u x ) will be given in values normalized to 
the unit cell length (u.c). As the nanowires considered here present several polymorphs, 
when ambiguous the reciprocal coordinates will be given a ZB (zincblende) or WZ 
(wurtzite) subscript, e.g. {Ill) zb and (002)wz which correspond to the same reflection 
with different reference lattices. 

2. Coherent Bragg Imaging on strained nano-objects: theory 

In this section, we present the general theory of Bragg coherent X-ray diffraction applied 
to strained nano-objects, and particularly the different approximations that can be used 
to compute the scattering, depending on the type of structural model that is available. 

In the case of focused X-rays, the wavefield close to the focal point is not a strict 
plane-wave, due to the curvature of the incident wave which induces a variation in both 
the amplitude and the phase of the X-ray wavefield. However, if we assume that the 
scattering object is close enough to the focal point and is small compared to the size of 
the focal point, then a plane wave is a quite good approximation [ID] . 

For small objects (<200 nm) it is reasonable to evaluate the scattering within the 
kinematical approximation. However, in the case of larger objects with heavy materials, 
refraction effects must be taken into account [20] . They will be ignored in this article 
(focused on nanowires with a diameter smaller than 100 nm) for the sake of simplicity. 

Given these approximations the scattering of the X-rays is due to the interaction 
with the sample electrons, and can be written as the Fourier transform of the electronic 
density: 

A(k)= f p(r)e 2 ^ = FT [p(r)] (1) 

where A is the complex scattered amplitude for the scattering vector k = kf — hi, p(r) 
is the electronic density in the sample, and FT denotes the Fourier transform. 

This formula is suitable for the analysis of CDI experiments in the (relatively) 
small angle regime, i.e. for small k values: indeed the real-space resolution at which 
the electronic density must be described is directly related to the extent of the Fourier 
transform in reciprocal space. For example in the case of a nanocube of silicon with a 
size of 100 unit cells (54.3 nm), in order to compute the scattering up to a (moderate) 
0.1 nm resolution, the size of the electronic density array would be 543 3 ~ 160 x 10 6 
points. 

This approach is unsuitable for the study of strained nano-objects using CBI, 
where the scattering will be collected far from the small angle regime, especially as 
the displacement from a perfect periodic lattice will often be a fraction of the unit cell's 
dimensions. In this case it is preferable to use an atomistic description of the crystal to 
compute its scattering, i.e.: 

A{k)=Y J U{k)e 2 ^ (2) 
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Figure 1. Coherent diffraction imaging of a single Si nanowire: (a) simulated 2D 
scattering pattern (in the plane perpendicular to the nanowire axis), (b) experimental 
pattern recorded on a 95 nm silicon nanowire |23j , and (c) the corresponding real-space 
reconstruction of the nanowire cross-section 



where fi(k) and fl are respectively the scattering factor and the position of atom i. In 
the case of a strained structure, the position of atom i can be written as: 

fi = rf + Ui (3) 

where r° is the ideal position of atom i in the unstrained crystal and Ui its displacement 
from that position. If the crystal is large compared to the size of the unit cell, and if 
the deformation varies slowly (i.e. all elements of the strain tensor are <C 1 anywhere in 
the crystal), the displacement in the crystal can be described using block unit cells (jj|): 

rXj = R°j + Uj + fl (4) 

where ry denotes the absolute position of atom i in unit cell j, fl is the position of the 
atom i relative to the unit cell, and Uj are the ideal position and the displacement 
vector of unit cell j. The scattering can then be rewritten as: 

A(k) = E E /i® e 2i7r ^+^) = F(k) £ e 2i7r ^4+^) (5) 

i j j 

where F(k) = Y,i fi(k)e 2t7Tk ' n is the structure factor of the crystal, i.e. this is only valid 
if all unit cells present the same crystallographic components. If ||(A; — k ) • Uj\\ 1 
(where ko is the scattering vector at the position of the Bragg reflection) then we can 
approximate k ■ Uj by k ■ Uj and we obtain the Fourier transform: 

A(k) = F(k) e 2l ^e 2m ^ w F{k)FT [ e 2i ^] (6) 
j 

This last approximation is generally used for the analysis of CBI experiments 
[H 1211 |22l |3T] as it allows a fast Fourier transform computation. However, its validity 

J This is only possible if the calculation is made around a symmetry-authorized Bragg reflection. In 
the case of a forbidden reflection, small relative atomic displacements within the same unit cell, as well 
as partially filled unit cells, can induce strong variations of the intensity of the 'forbidden' reflection 
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strongly depends on the size of the object and the amplitude of its deformation: e.g. 
for an object of size 100 unit cells, the extent in reciprocal space of its scattering will 
be N * Y5o> w here iV is the number of oscillations recorded. For TV = 10, and if the 
maximum deformation on the sides of the crystal is equal to one unit cell (corresponding 
to a 1% strain over 100 unit cells) then the maximum difference between k ■ Uj and ko ■ Uj 
will be equal to 0.1, which cannot be neglected (2% ■ 0.1 ~ 0.63 rad w 35°). In such 
a case it is more accurate to use directly equation (j2J) (if individual atomic positions 
are known) or equation <^ (if the displacement field is described per block unit cell), 
although the calculations are significantly slower. 

Finally, it is important to note that the strain sensitivity of a single CBI 
measurements is limited to the direction parallel to the scattering vector: a successful 
object reconstruction will yield 2nk- U as the phase of the reconstructed density. A full 
reconstruction of an object deformation will require data collection around three linearly 
independent reflection scattering vectors for the same object [251 E2], less if symmetry 
can be used. 

3. Coherent Bragg imaging on single nanowires in the absence of faults 

3.1. Homogeneous nanowires 

Homogeneous nanowires without any external stress only present strain on the atomic 
layers that are closest to the surface [33], so given the limited real-space resolution 
obtained by CBI (usually at least 5 nm), this outer layer contraction should be negligible. 
In that case coherent scattering around a Bragg reflection will yield the same signal as 
for the small angle range, i.e. the Fourier transform of the shape of the nanowire. As 
the nanowire length is much larger than its diameter, the signal will then be a quasi- 
2D signal which corresponds to the 2D FT of the nanowire cross-section. Example of 
simulated and experimental CBI signal recorded on a single silicon nanowire is presented 
in figure [JJ (see [23J for details). 

As in all X-ray experiments, the phase of the scattered amplitude is lost during 
the experiment, and must be recovered either using ab initio phase retrieval algorithms 
[TH [TBI [TBI [T7] or using a model with some a priori (e.g. symmetry) information [23] . 

3.2. Heterogeneous nanowires 

In the case of an heterogeneous nanowire, CBI will be sensitive to the displacement and 
the nature of all atoms. We have simulated an [OOljiy^-oriented InAs nanowire, with an 
InP insertion (see figure E]), using atomistic simulations - both sections have a wurtzite 
structure. The radius of the regular hexagon is R=60 nm, and several wire lengths and 
insertion thicknesses were simulated (see figure |2]) to test the sensitivity of CBI. The 
relaxed atomic positions are computed using Keating's valence force field (VFF model) 
[Ml |35]. 
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Figure 2. Atomistic simulations of an InAs wurtzite nanowire (0 = 60 nm, height= 
100 nm) with (a,b,c,d) a 3 nm InP insertion and (e,f,g,h) a 15 nm insertion: (a,e) the 
radial displacement (along [110], expressed relatively to the perfect InAs lattice, in 
unit cells (u.c.)), (b,f) the axial displacement (along [001]), the calculated intensities 
around the (c,g) (004) and (d,h) (202) reflections, with a logarithmic colour scale. The 
map coordinates are given in reciprocal lattice units (r.l.u.) relative to InAs. 

The amplitude of the displacement field (computed relatively to a perfect InAs 
lattice) has a maximum value of 1-2 unit cells in the horizontal direction (u x ), depending 
on the thickness of the insertion; in the vertical direction (u z ) the maximum displacement 
varies with the InP thickness t, ranging from 0.16 (£=1.5 nm) to 0.6 unit cells (£=15 
nm). The InP insertion has a lower bulk lattice parameter than InAs (dyfiz = 0.415 nm, 
a wz s = 0.4268 nm), so that the insertion contracts the lattice, as can be seen in figure 

m 

Due to the smaller lattice parameter, scattering from the InP insertion occurs at 
higher angles: although a clear perturbation in the scattering with t=3 nm can be seen 
(figure [2Tc) and (d)) - an homogeneous nanowire would present roughly (jjj) the same 
signal around the two reflections-, it is only with a larger insertion (figure [2(g) and 
(h)) that a well-separated peak appears for the InP insertion, roughly two orders of 
magnitude less intense than the main InAs peak. 

§ The only difference being due to the amplitude of the atomic scattering factors 
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These simulations indicate that it is possible to measure a scattering signal specific 
to a small (a few nm) insertion, assuming at least 3 or 4 orders of magnitude in the 
experimental data, which is usually achieved. However, this simulation assumes a perfect 
nanowire, without any structural fault, which is difficult to achieve. It is therefore 
important to evaluate the influence of faults on the coherent scattering, and whether it 
can hide the scattering signal due to a strained lattice. 

4. Coherent Bragg Imaging on nanowires with stacking faults 

4-1- Stacking faults in nanowires 

Stacking faults commonly occur in nanowires [361 E3 EEJ ESI HQj EJ H2, 03] and some- 
times lead to a periodic twinning structure [HI HSl HH1 STJ HH] . This is especially true 
for III-V nanowires, where the nanowire geometry and growth kinetics may stabilize the 
metastable wurtzite structure instead of the zincblende bulk phase [49j [501 EE] • 

This section presents first a theoretical study on the influence of stacking faults 
when studying a faulted nanowire, and then shows how it is experimentally possible to 
avoid the effect of the stacking faults when studying strain in a nanowire. Finally we 
show how it is possible to use CBI in order to retrieve statistical information about a 
fault sequence in a nanowire. 

4-2. CBI in the presence of stacking faults: simulations with InAs/InP nanowires 

During a coherent diffraction experiment, the scattering around a Bragg peak is sensitive 
to any displacement from a perfect periodic lattice. For example in the case of a 
wurtzite nanowire, a simple deformation fault [29J corresponds to the following sequence: 
AB|AB|AB|AB|CA|CA|CA|CA|CA, with the first part being the normal 'AB' sequence, 
followed by an 'AC sequence, where all atoms atoms are shifted by v = |(110)vkz 
from their position in the normal sequence. This abrupt shift will produce a specific 
interference pattern around reflections for which the scattering vector is not orthogonal 
to v, and these interferences will cumulate with those coming from the strain-induced 
deformation. 

An example of simulated scattering from InAs/InP nanowires with one and three 
stacking faults is given in figure [3j in both cases the oscillation fringes are affected by 
the faults. In theory, such faults would not prevent the reconstruction of the density 
and displacement fields in the wire by ab initio phase reconstruction algorithms, as the 
faults correspond to a simple 'phase shift' of e 2twk ' u (e.g. for a k = (202) reflection and 
v — |(o— b), the phase shift is 2irk ■ v = ^-). 

However in practice, the reconstruction would be hindered by the fact that the 
resolution of the experiments is currently limited to 5-20 nm (e.g. ~ 10-40 unit cells 
in the case of InAs), and therefore a displacement configuration such as seen in figure 
El^b) may be missed. Moreover, reconstruction algorithms for strained objects tend to 
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Figure 3. Simulation of the influence of stacking faults on the coherent scattering 
from InAs/InP nanowires (0 = 60 nm, height= 100 nm) around a (202) wz reflection. 
2u x + 2u z for (a) the original InAs/InP nanowire obtained using atomistic simulations, 
and in the case of (d) one and (e) three growth faults. The color scale is expressed 
in unit cells (u.c.) and is the same for (a,d,e). The resulting scattered intensity 
is shown in (b) (one fault) and (c) (three faults); these can be compared to figure 
HJh) (no fault). The map coordinates are given in reciprocal lattice units (r.l.u.) 
relative to InAs. There are a number of small changes from the scattering from 
the original InAs/InP nanowire (see figure [2jh)), notably with shifted fringes in the 
domain H e[2.00; 2.04], L e[1.95, 3.10]. The effect of the faults on the reconstruction is 
more visible when computing the inverse Fourier transform of the simulated complex 
amplitude around the (202) reflection: the complex recovered density (see text for 
details) is shown in the case of one (f) and three (g) faults. In these images the phase 
{2irk ■ U) is given by the colour (as indicated by the colour wheel), while the amplitude 
is given by the saturation of the colour (the saturation is given on a linear scale: 
white corresponds to no electronic density and full saturation to the normal average 
zincblende density). The reconstruction yields in both cases a low density (white) at 
the position of the structural faults. The inserts show a plot of the amplitude along 
the vertical axis in the center of the reconstructed wire, where the dips in intensity 
can be clearly seen. 
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converge more slowly than in the unstrained the real-space constraint in the 

reconstruction can only use the finite support for the object and not its positivity ; for 
this reason it can be useful to impose continuity to the phase field of the reconstructed 
object [52], a strategy that should be modified to take stacking faults into account. 

Moreover, as already pointed out in [22], the scattering around a given reflection 
will depend on the domains (zincblende or wurtzite) that are found along the wire: 
each type of domain may only contribute to the scattering of reflections allowed by its 
structure, as will be detailed in section 14.41 

Finally to quantitatively evaluate the influence of the stacking faults of the nanowire 
reconstruction from CBI, it is interesting to use the simulated data (figures [3] (b) and (c)) 
and compute the inverse Fourier transform ([jj]) to obtain pe 2l7T ^'^ . In order to correctly 
recover one (xOz) plane in direct space, this was performed from 3D simulated data 
around the (202) wz reflections presented in figure [3j The result is presented in figure 
[3]T) and (g), in the case of one (f) and three (g) inserted growth faults. The deformation 
field is recovered as the phase of the signal (corresponding to 2ir(2u x + 2u z )), but the 
most striking features are the 'holes' in the reconstructed amplitude at the position of 
the faults. 

These holes are an artefact caused by the presence of the faults, and can be 
simply explained: they introduce an (ABC) succession of layers in the middle of the 
(AB|AB|AB) sequence, and this particular sequence leads to a null intensity. If the first 
layer (A) is used as a reference, the second (B) is shifted by vb — v + (00|) = {\\\)wz 
and the third by vq = 2v + (001) = (||1)vkz- Their scattered amplitude around the 
(202) wz reflection is then equal to: 

1 + e 2M B + e 2ink-v c = 1 + e if- + g Sf = q ^ 

As a consequence, this succession of layers of the nanowire contribute no intensity 
in this part of the reciprocal space, and it logically follows that the inverse Fourier 
transform from the scattered amplitude around the (202) reflection yields a null 
electronic density near the faulted parts of the wire. 

Recovering the correct intensity near the faults would require a direct space 
resolution sufficient to resolve individual layers, i.e. half a wurtzite unit cell along 
the (001) direction. This corresponds to a range in reciprocal space such that AL > 2, 
i.e. at least including the (201) and (203) reflections. In practice, it will be much easier 
to select a reflection such that k ■ v is integer and therefore will be insensitive to the 
faults. An example of different reconstructions using several reflections from a faulted 
structure can be seen in [32] . 

|| The following method was used: the direct (real space to reciprocal) calculation used was performed 
using atomic positions as denoted in equation [21 The points in reciprocal space at which the 
scattering was computed were located at the following reciprocal lattice unit coordinates: h e[1.76, 2.4[, 
k e[— .24, .24[, I e[1.76,2.4[, with a step equal to .004. The reconstructions (reciprocal to real space) 
where computed using a Fast Fourier Transform, and the resulting array was cropped to keep only the 
area containing the crystal - only a 2D section is presented in figure El 
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Figure 4. (a) Transmission electron microscopy image of an InSb/InP interface, 
with a twinning fault in the bottom of the InSb part. Experimental diffraction image 
((333) zb reflection) from an InSb/InP nanowire, measured at the (b) top and (c) 
bottom (near the InSb/InP interface) of the InSb section. The lattice parameter of 
the InP section is sufficiently different (10.3%) not to diffract on the same image, (d) 
Simulated deformation as a function of the height near the base of the InSb section 
(height=250 nm), and (e) the corresponding scattering simulation, which matches well 
the measurement recorded at the bottom of the wire. 



4-3. CBI in the presence of stacking faults: measurement of one InSb/InP nanowire 

In order to measure the strain from nanowires in the presence of known stacking faults, it 
is also possible to select a reflection which is insensitive to the displacements associated 
to the faults, e.g. for a displacement vector v, any Bragg reflection for which k ■ v is an 
integer will be unaffected by the fault. 

We tested this approach on InSb/InP heterogeneous nanowires. These nanowires 
were grown [531 l5~4"l 155] from Au seed particles with 30 nm diameter on InP(lll)s 
substrates using metal-organic vapour phase epitaxy. The InP segment was first grown 
at 420°C after which the InSb segment was grown at 450°C for 13 min using a V/III 
ratio of 70. As can be seen in figure IH(a), in these wires the InP section has a wurtzite 
structure with many stacking faults - this part was not studied to avoid diffraction from 
the InP substrate. The InSb section has a zincblende structure, with a few (one to 
three) twinning faults (j3k) located near the base, were there is also a strain relaxation 
in the first 10-20 nm. 

We measured the scattering from single nanowires on the ID01 beamline of 
the ESRF, using a Fresnel-Zone-Plate (FZP) to focus the 8 keV beam down to 
300 (vert.) x 500 (horiz.) nm 2 (FWHM), which allowed to select different sections 
of the nanowire for diffraction (the length of the InSb section is 0.68 ± 0.15 fim). As 
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shown in figure H(b) and (c) we measured coherent diffraction from both the top and the 
bottom of the InSb (the sample was vertically translated by 250 nm between the two 
images). The image corresponding to the bottom of InSb features additional fringes, 
which cannot be due to twin faults as the chosen reflection -(333) ZB - is insensitive to 
stacking faults. Note that the 'top' image (figure 11(b)) also features an asymmetric 
shape with some diffuse scattering at larger A; <m> values, but that is probably due to 
some contribution from the base of the wire due to the tails of the focused X-ray beam. 

In order to reproduce these asymmetric fringes, we simulated the bottom of the InSb 
nanowire (height=250 nm, width=54 nm), where the spacing between successive layers 
(see figure 11(d)) is increased by e° zz = 2% at the bottom, with a simple exponential 
decrease e zz the precise relaxation law could not be predicted (the 

relaxation due to the 10.3% lattice mismatch occurs mainly through dislocations and 
radial relaxation), we tried several values for A, with A = 10 nm giving the best 
agreement with the experimental pattern. 

The simulated scattering (figure 11(e)) reproduces well the fringes, which are due to 
the combination of the height of the wire selected by the beam (the fringes period) and 
the deformation (the asymmetric intensity). The shorter dm-spacing at the bottom of 
InSb (relatively to bulk InSb) is probably due to interdiffusion during the deposition, 
with an InSb x Pi- x chemical composition varying in the first 10-20 nm of InSb (see 
figure 4(c) in [55]). 

While this experiment has shown that it is possible to be sensitive to strain near 
an heterogeneous interface (even in the presence of faults), a new measurement with 
more intensity and extent in reciprocal space would be necessary to enable the ah initio 
reconstruction of the deformation. 

4-4- CBI with many stacking faults: measurement of a GaAs/GaP nanowire 

4-4- 1- Experimental measurement: In the case of a nanowire with many stacking faults, 
Bragg coherent diffraction can be used to get information about the fault sequence. We 
studied GaAs/GaP nanowires on beamline 34ID-C of the Advanced Photon Source; 
during this experiment we used a 1.1 x 2.2 /im 2 beam, obtained with Kirkpatrick-Baez 
mirrors [7j. Slits were closed down to 20 x 50 [im 2 before the KB mirrors to select 
the coherent part of the beam. The scattered intensity was collected with a direct 
illumination CCD camera with 22 /j,m pixel size, placed at 1300 mm from the sample. 

The samples (see [38] for synthesis details) consisted of nanowires divided in two 
parts, with GaP on the bottom and GaAs at the top, and an average diameter of 50 
nm (figure [5(a)). They were grown on Si02 and were randomly oriented. In order to 
collect the diffraction from the GaAs section (presenting fewer faults), the detector was 
placed at the 29 angle corresponding to the (111) reflection, and then the samples were 
scanned until a diffraction pattern appeared on the detector. 

For some wires the recorded images did not present the expected Fourier transform 
of the nanowire cross-section (figure [1(a)), but rather presented a 'bar-code' pattern, as 
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presented in figure E(b): from this image the width of the nanowire can be estimated 
to 50 nm from the horizontal width of the bar-code pattern, while the width along the 
length (vertical) of the barcode pattern indicates that it is related to domains along the 
wire which are about 10 nm long. 

These measurements performed on a single NW can be compared to those obtained 
using non-coherent X-ray scattering on an assembly of epitaxial wires: as an example, 
figure [5](e) shows the experimental intensity recorded along (20L) WZ crystal truncation 
rods (CTR), on the same InSb/InP nanowires as described in section 14.31 and figure 
Ufa). The grazing- incidence diffraction was recorded on the BM32 beamline of the 
ESRF, with a 0.05° incident angle to minimize the contribution from the InP substrate. 
The two CTR are extremely different, due to the fact that the InSb part (black 
line) has a zincblende structure with few faults (jjj) (e.g. the sharp peak around 
(20|)vkz = (311)zb)? while the InP section (red curve) is predominently wurtzite with 
many defects, and only features much larger peaks. Note that this type of CTR recorded 
on an assembly of nanowires can be analysed using statistical models of the stacking 
sequence [29] . 

4 -4 -2. Model for the faulted nanowires: In the case of < 111 > nanowires with 
stacking faults, and for the (111) zb = (10|)vkz reflection (j+|) , the phase shift of the 
scattered signal is e 2lnk ' v , where vwz = =t|(^ — Depending on the shift from an 
arbitrary reference stacking, each domain along the wire then diffracts with a phase 
shift if G [— yjO,+yL i- e - the nanowire is seen as a 'phase object', and the scattered 
signal is equal to the Fourier Transform of e l,p . 

The configuration along the wire can be relatively complex: while GaAs as a bulk 
crystal is stable in the zincblende structure, GaAs nanowires can either be wurtzite 
[4"9] or zincblende (see figure 2(c) in [38]). They can therefore be expected to have 
4 types of domains, with two possible twins for zincblende and wurtzite, and each of 
these domains can be shifted by v or 2v from a reference configuration. The scattering 
from these different domains is however relatively simple, as only one zincblende twin 
configuration will contribute to the scattering: the (ill) zb reflection corresponds to the 
(10|) in wurtzite reciprocal lattice coordinates, and therefore has no intensity. As for 
the twin zincblende structure (which is the image of the normal zincblende with respect 
to a mirror perpendicular to the < 111 > direction), the corresponding reciprocal lattice 
coordinates are |(511), which does not exist either. As a consequence the reconstructed 
domain structure from the CDI data should yield both the type of domain, diffracting 
(normal zincblende) or not (twin zincblende or wurtzite) and the corresponding phase 
shift of these domains. 

% There are two twin variants of the zincblende InSb, but as one wire is of a given variant (with up to 
to three twin faults), it only contributes to one CTR and therefore does not enlarge the Bragg peak. 
See figure 3 in [35] for a description of the scattering of the different variants. 
+ Using as basis vectors in direct space: awz = ^(llO)zs, t>wz = f (OH)zb and cwz — (H1)zb 
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4- 4-3. Domain reconstruction: In order to reconstruct the phase sequence from the 
scattering data, it is theoretically possible to retrieve the phase from the scattered 
signal using phase retrieval algorithms (and then using an inverse Fourier Transform 
to reconstruct the object), but this requires that the data be sufficiently oversampled. 
In our case, the barcode pattern was very sharp (see figures Mjo) and (d)), with some 
peaks as small as 2 pixels, and the wire was longer than the beam. As the number 
of parameters was limited (400 integrated data points), we decided to retrieve the 
phases of each domain in direct space, using parallel tempering, a biased reverse Monte- 
Carlo algorithm [56]. During each cycle of this algorithm, either a small number of 
phases or an occupancy factor (1 or 0) along the wire are randomly changed. The new 

configuration is then accepted if the R-factor (R = J V^/p °\2 C ) diminishes, and 

otherwise the configuration is rejected with a probability P = e~( R ~ Rprevious ^ T , where T 
is the temperature of the algorithm. 

The object itself was modelled as a ID object with rj 5 nm domains, each with a 
phase ip G [— 0, +y) an d with an occupancy (1 for a diffracting domain, otherwise). 
As the GaAs part of the wire was longer than the beam at the focal point, the beam 
was modelled with a gaussian shape, with a 1.1 [im beam size. The reconstructed phase 
object is presented in figure E^c). Note that the reconstructed object is not unique, 
which is a known issue in the case of ID problems [57]. However the analysis of several 
solutions show that they have the same statistical information (number and size of each 
type of domain), as the uncertainties given below indicate. 

Quantitative information can be obtained from the reconstruction. First, the size 
distribution of the normal zincblende domains (the only ones seen by the reconstruction) 
is as follows: 5 nm (54 ±1%), 10 nm (28 ± 1%), 15 nm (8 ± 4%), 20 nm (8 ± 4%), with 
an average domain size of 9 ± 0.5 nm. This value is comparable to the size that can be 
estimated from the FWHM of the barcode pattern (figure E(d)), L = — 11 nm. This 
length is also in agreement with the apparent frequency of 'kinks' on figure EI a), which 
are known [451 H6j HHJ H2] to originate from a succession of zincblende twins. 

Furthermore the results of several (10) optimization runs indicate a per- 
centage of occupied (normal zincblende) domains slightly lower than 50%: 
P (normal zincblende) = 46 ± 3%. As the two zincblende twins should be equiprobable 
(any bias introduced by the terminating planes of the GaP section should vanish after a 
few domain changes), this indicates that the GaAs section consists of 92 ±6% zincblende 
domains, and therefore of 8 ± 6% wurtzite domains. 

4-4-4- Discussion on improvement to the ab initio domain reconstruction: Structures 
like the 'bar code' pattern are found to be particularly difficult to phase and this has not 
yet been possible with traditional methods [T3J, [15] . There are several possible reasons 
for this: first and foremost is that the interference pattern is only one-dimensional (ID), 
even though the full recorded pattern is 2D. As mentioned above, the cross-section 
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Figure 5. (a) Scanning electron microscopy image of a < 111 >zb nanowire with 
a GaP bottom section and a top GaAs section. Note the nano-facetting of the wire, 
which can be related to stacking faults along the length of the wire, (b) Coherent 
diffraction image recorded on a (lll)zs GaAs reflection. The "barcode" pattern is 
due to the faults: the nanowire can be modelled as a stacking of domains with ±^ 
phase shifts (see text for details), (c) Reconstructed phase object from the diffraction 
data: each domain along the wire corresponds to a phase <p 6 [— 0, ]- associated 
with colors red, green and blue, or white for a domain not contributing to the scattering 
of the (111) zb reflection. The intensity of the color corresponds to the amplitude of 
the beam (gaussian with FWHM=1.1 fim) along the wire, (d) Integrated ID profile 
of the 'barcode' pattern, observed (red) and simulated (black) using the phase object 
in (c). (e) Non-coherent crystal truncation rods (CTR) measured on an assembly of 
InSb/InP nanowires: the CTR from InSb (black), which has a zincblende structure 
with very few faults, presents very sharp peaks, whereas the InP CTR (red) is very 
wide, due to the high density of stacking faults in the wurtzite structure. 



is just the shape transform of the nanowire itself; without the faults the diffraction 
pattern would be a single peak, being a disk viewed edge-on. ID diffraction patterns 
are prone to factorisation of their amplitude functions, so the object can be expressed 
as a convolution of substructures; 2D and 3D functions are much less likely to factorise 
[57] . The complex conjugate of each of these pieces gives the same modulus for its 
factor, so there are a large number of equivalent solutions [57] . 

The second problem, somewhat related to the first, is that the sequence of nanowire 
faulted segments tends to be self similar, with internal repeats; permutations of these 
repeating subsequences give virtually the same diffraction pattern, so the solution is 
non-unique, especially in the presence of noise. This non-uniqueness was indeed seen 
with the biased Monte Carlo method [56J described above, just as was also found with 
the 2D FeAl antiphase domain inversion results reported previously [19J. 

A possible solution to the uniqueness problem has been proposed by Rodenburg 
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Figure 6. Ptychography data obtained by moving the focussed beam in 1-micron 
steps, recorded for the same GaAs/GaP nanowire as shown in figure [5j The sequence 
of diffraction patterns shows a steady evolution of the positions of the maxima and 
minima. 



et. al. in a method called "Ptychography", in which a well-defined beam illumination 
function is scanned over the region of interest in the sample, along the wire in this case. 
The overlapping diffraction patterns that result contain redundant information and the 
overlap can be made arbitrarily small. An example of ptychographical data, recorded 
for the same GaAs/GaP nanowire, is shown in figure |6] ; the individual diffraction 
patterns are clearly related to each other, but show differences that contain additional 
information about the structure. An algorithm for the recovery of the structure has 
been proposed [58J, but not yet tried with these data. 

Partial coherence is another important issue that may limit the effectiveness of the 
phasing methods described here. Typically, the size of the objects discussed here is 
considerably smaller than the lateral coherence length of the X-ray beams used. Yet 
partial coherence effects can be seen directly in the graph of figure Efd): the visibility 
of the data (red) is less than that of the best-fit simulation (black). The data rarely 
come down near zero intensity, while the simulation comes closer; this is because there 
is a second intensity component present in the data, which does not interfere with the 
main one. Given that the limited lateral coherence originates from the finite size of the 
synchrotron X-ray source, it can be modelled as a Gaussian smearing of the resulting 
diffraction patterns [59J, but this does not lend itself to easy methods of removing it. 

A more interesting way forward is offered by the new work of Flewett et. al. [60] . 
who show that a typical synchrotron source can be modelled by a small but finite number 
of modes, with more than 90% in the primary mode in the case studied. The additional 
modes can be attributed to the edges of the source, particularly in the horizontal di- 
rection. These extra modes, with an illumination profile that is orthogonal to the main 
mode, give different diffraction patterns from the same object; these are combined (as 
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intensities) in the diffraction pattern which therefore has lower contrast (visibility). So 
long as the extra modes are a small fraction of the main mode, they can be used to 
correct for the distortions of the data within an algorithm, without adversely affecting 
the iterative computation on the main mode. Again, this method has not yet been tried 
on the data of figure El 

5. Conclusion 

During the last ten years, coherent X-ray scattering has seen very significant progress, 
notably through the use of better focusing optics, producing a smaller beam size with 
a much higher photon flux (5.10 3 — 5.10 4 ph/s/nm 2 ). This allowed to study objects 
smaller than 100 nm, both for small-angle and Bragg CDI experiments. 

The development of Bragg CDI is particularly important as it gives access to 
the deformation in the volume of objects at the nanoscale, which is essential for the 
understanding of their structural properties. This point has been illustrated by our 
experimental and simulation studies on InAs/InP and InSb/InP nanowires. In the case 
of faulted structures, CBI can be used to retrieve a statistical image (size and type of 
domains) of the wire. The most important development is probably the ability to use 
the very small size of the beam to study different parts of a given heterostructure. 

However as we have shown in this study, information about the shape and strain 
state can be hindered by the presence of structural faults, such as stacking faults in 
nanowires - as we have shown this can lead to voids in the reconstructed electronic 
density. While CBI allows a quantitative analysis of the structural faults, illustrated 
here in the GaAs system, the combined presence of strain and stacking faults in a 
structure remains a complicated problem - in the simplest case (stacking faults) the use 
of carefully selected reflections will allow to be insensitive to the stacking order, but 
this requires the ability to study several reflections for the same object, which is still 
difficult for objects smaller than 100 nm, as the mechanical precision (confusion sphere) 
of existing goniometers is much larger than this value, especially when performing several 
rotations to select different reflections. 
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